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Abstract. We study the asymptotic behavior of the heat content on a com- 
pact Riemannian manifold with boundary and with singular specific heat and 
singular initial temperature distributions imposing Robin boundary conditions. 
Assuming the existence of a complete asymptotic series we determine the first 
three terms in that series. In addition to the general setting, the interval is 
studied in detail as are recursion relations among the coefficients and the re- 
lationship between the Dirichlet and Robin settings. 
Subject classification: 58J32; 58J35; 35K20 



1. Introduction 

1.1. Historical framework. Let M be a compact m-dimensional manifold with 
boundary DM and Riemannian metric g. Let D be an operator of Laplace type 
which drives the evolution process (see Section IT~2l for details). Let <fi represent the 
initial temperature of M and let T{x; t) be subsequent temperature of the manifold. 
A good general reference to the heat equation is provided by [16j HB] ■ We impose 
suitable boundary conditions B to ensure the process is well defined. Let p be the 
specific heat. We suppose for the moment </> and p are smooth and let /3(t) be 
the heat content of the manifold. There is a complete asymptotic series for j3 as 
t I with locally computable coefficients. We postpone for the moment a precise 
definition to avoid complicating the exposition at this stage. 

The problem was first studied by [2J [3] in the context of domains in Euclidean 
space with smooth boundaries. We note that regions with polygonal boundaries 
were examined in [5]. Subsequent results for regions with fractal boundaries were 
obtained in [TJ [3J. Apart from these papers, most of the work on the heat con- 
tent asymptotics has been in the smooth category and we shall always assume the 
boundary of M to be smooth. The "functorial method" has been used to express the 
asymptotic coefficients in terms of geometrical quantities for Dirichlet and Robin 
boundary conditions by [5J El HE UU ■ Other authors [H EH EE 121 130] have ex- 
amined these coefficients using other methods. Various other boundary conditions 
also have been considered [13j [2Gj [21] as have variable geometries [19] . Motivated 
by work of [15] for the growth of the heat trace asymptotics, the growth of the heat 
content asymptotics has been examined Ell [551 GUJ ■ 

The case where the initial temperature and the specific heat are singular in a 
controlled fashion at the boundary will form the centerpiece of this paper - we 
refer to Section 11.31 for precise definitions. A rigorous treatment was provided in 
[14] where the initial temperature is singular on the boundary but the specific heat 
is smooth. Subsequently, the case where the specific heat can be singular as well 
was studied for Dirichlet boundary conditions |1Q[ 113) . The present paper examines 
the situation for Neumann boundary conditions. 

1.2. Operators of Laplace type. Let (M, g) be a compact Riemannian manifold 
of dimension m with smooth boundary DM. Let V be a smooth vector bundle over 

l 



2 



M. VAN DEN BERG, P. GILKEY, AND H. KANG 



M and let D : C°°(V) — > C°°(V) be a second order partial differential operator 
on the space of smooth sections to V. We adopt the Einstein convention and sum 
over repeated indices. We say that D is of Laplace type if the leading symbol of D 
is given by the metric tensor, i.e. locally D has the form: 

D = -{g^d x ^ Xv +A°d Xa +B) . 

The following Bochner formalism |18j permits us to work tensorially: 

Lemma 1.1. There exists a unique connection V onV and a unique endomorphism 
E ofV so that D = D{g,V,E) = - (g^ Va x „ Vg +E). The connection 1-form uj 
of V and the endomorphism E are given by 

w P = \{9nvA v + g ae T aEt , Id) andE = B- g^ v (d Xu w„ + Ull w v - uj a T^ a ) . 

If D — d*d + dd* is the Laplacian on the space of smooth p- forms, then V is 
the Levi-Civita connection and E is given in terms of curvature. If D is the spin 
Laplacian, then V is the spinor connection and E is a multiple of the scalar cur- 
vature. We shall use the Levi-Civita connection and V to covariantly differentiate 
tensors of all types. 

We must impose suitable boundary conditions. Let e m be the inward unit normal 
vector field on the boundary. Let S-r be an auxiliary endomorphism of y|aM and 
let V be the connection on V given in Lemma 1 1.1 1 Let Bx> and Bn be the Dirichlet 
and the Robin boundary operators which are defined, respectively, by setting: 

Bv4> = 4>\dM and Bn^={V em 4> + Sn^)\dM for feC°°{V). 

The Neumann boundary operator is defined by taking S-jz — 0. We let Dp / K denote 
either the Dirichlet or the Robin realization of D henceforth. It is convenient to 
have a common notation. 

1.3. The initial temperature and the specific heat. Let r be the geodesic 
distance to the boundary. This is a continuous function on all of M which is 
smooth near dM. If y = (y 1 , ?/ m_1 ) is a system of local coordinates near a point 
P € dM, then x = (y,r) for r € [0,e] is an adapted system of local coordinates 
near P in M for some e > 0; the curves 7 y (s) := [y, s) are unit speed geodesies 
perpendicular to the boundary and d r is the inward pointing unit geodesic normal. 

To have a common notation, let W = V or let W = V*. If a € C, let K a (W) be 
the vector space of all sections $ to W which are smooth on the interior of M such 
that r Q $ is smooth near dM. The parameter a controls the blow up (resp. decay) 
of $ near the boundary if 3t(a) > (resp. SR(a) < 0). Let (j) £ JC ai (V) represent 
the initial temperature of M and let p £ JC a2 (V*) give the specific heat. We shall 
always assume that 

9£(ai) < 1, 3?(a 2 ) < L and at\ + a 2 % . (l.a) 

The first two conditions ensure that <f> and p are in L . The final condition is 
imposed to avoid interactions between the interior and boundary terms in the as- 
ymptotic series for the heat content. We will return to this point subsequently. 
Expand <j> £ (V) near dM in a modified Taylor series: 

oo 

cf>{y,r) ~r- ai £y(y)r>' as r ± 0+ where <fP = jri^o^ {r ai 4>)\dM ■ (Lb) 

3=0 

We shall usually be working with scalar operators and can choose a local section s 
so that V em s = 0. We may then regard <j> as a function and the above expansion 
as a Taylor series. A similar expansion, of course, is valid for the specific heat p 
where we regard p as a section to the dual bundle and covariantly differentiate with 
respect to the dual connection on V* which has connection 1-form —uj*. 
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1.4. The heat equation and the heat content. If <fr 6 K. ax (V), then the sub- 
sequent temperature T := e~ tDB 4> is characterized by the relations: 

(dt + D)T = (evolution equation), 

limt_|,o T(-; t) — cf> (initial condition), (Lc) 

BT(-;t) = for t > (boundary condition). 

Let (•, •) be the natural pairing between V and the dual bundle V* , let dx be the 
Riemannian measure on M, let dy be the Riemannian measure on dM, and let 
p € /C Q2 (V^*) be the specific heat of the manifold. The total heat content of the 
manifold is defined by setting: 

P(4>,p,D,B){t) := / {e- tD ^,p)dx. 
Jm 

There is a smooth heat kernel K = Kd,b so that T(x;t) — J M K(x, x; t)4>{x)dx, 
and 

P(<t>, p, D, B)(t) = / (K(x, x; t)(f>(x), p{x))dxdx . 

JMxM 

This is well defined for <fi 6 ^(V) and p E L 1 (V^*); it was for this reason that we 
assumed 3fJ(or) < 1 and 3?(a 2 ) < 1 in Equation (jl.aj) . 

If is self-adjoint, as will be the case for either the Dirichlet or the Robin 
realization of the scalar Laplacian, then we can take a spectral resolution {0„, A n } 
for Db- Here the {</> n } is an orthonormal basis for L 2 (V) such that D<f) n = X n (f> n 
and B(f) n = 0. We then have 

K(x, x;t) = ^2 e^ tXn (f> n (x) <E) <p n (x) . 

This series converges in the C°° topology for t > 0. The temperature T is smooth 
for t > 0. However, the convergence to <f> as 1 1 in Equation (jl.cj) is not pointwise 
but in L l (V). 

The following example is instructive. Let M — [0,7r] and let D = —d%. The 
spectral resolution of the Dirichlet Laplacian is given by {n 2 , \/2/7r sm(nx)} n >i. 
Thus if (/) — 1, we have 

Tlx; t) = - V -e"*™ 2 sin(na;) . 
7r f— ' n 

n— l,n odd 

This is smooth for £ > 0. However, since T(0;t) = for t > and (j> — 1, the 
convergence is not pointwise. The associated heat content is given by: 

oo 

/3(1,1,A,^)(<) = - T ~ e ' tn2 
7r A — ' n z 

n— l,n odd 

= TT - ^t 1 / 2 + 0(t k ) as U V/c > 1 . 

1.5. The form of the asymptotic series. We shall need to consider certain 
integrals which are divergent and which need to be regularized. For example, the 
integral J M ((f>, p)dx is divergent if 1 < 5R(ai + 0:2) < 2. The Riemannian measure 
is not in general a product near the boundary. Since, however, dx — dydr on the 
boundary of M, we may decompose 

(<t>,p)dx = (c/) ,p )r- ai - a2 dydr + C^r 1 -" 1 -" 2 ) . 
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Let C £ := {x € M : r(x) < e} be a small collared neighborhood of the boundary. 
For + a 2 ) < 2 and a\ + a 2 7^ 1, define: 

l g Reg (<f>,p) := / {4>,p)dx+ f {(4 > ,p)dx~(cf > a ,p )r- a ^dydr} 
+ f (^V^xe^-^CL-ai-aa) -1 . 

JdM 

This is clearly independent of e and agrees with J M {4>, p)dx if 3t(ai + 0:2) < 1- The 
regularization lR g ((/), p) is a meromorphic function of cki + a 2 with a simple pole at 
u\ + a 2 = 1. More generally, the integrals (D n <p,p) which appear in Coniecture ll.2l 
below need regularization if 5R(ai + a 2 ) > 1 — 2n. Poles can appear whenever 
5R(o!i + a 2 ) = 1 — k. These poles are evident in the formulas given subsequently; 
we expect the local formulas for the asymptotic expansion of the heat content may 
involve log terms when a% + a 2 G Z and for that reason excluded these values in 
Equation (|l.a|) . 

To simplify the notation, we set: 

/3f ^,p,D) := (-lr/nl-l^UD^p)}. 

If M is a closed manifold, these are the invariants which would appear in the 
heat content expansion. We assume that the following conjecture (which extends 
the discussion of [10]) holds henceforth. We will justify the powers £( 1+J ~ Ql-Q2 )/ 2 
subsequently in Section 0] using dimensional analysis. We refer to [T3] where a 
related result was established when the specific heat is smooth. 

Conjecture 1.2. If (a\,a 2 ) satisfy Equation U.a\) . then there is a complete as- 
ymptotic series as t -\, of the form: 

00 

f3(cf>,p,D,B v/n )(t) ~ £>Af £>) 

00 

3=0 

The coefficients fif^ Q2 (4>, p, D, B-p /n) are given by integrals of local invariants over 
the boundary. 

Remark 1.3. This conjecture has been established in [T3] using the calculus of 
pseudo-differential operators in the special case that a 2 £ N or that ot\ € N. The 
extension to the present setting is motivated by that work. 

Let Rijki denote the Ricmann curvature tensor; with our sign convention, we 
have that -R1221 = +1 on the unit sphere S 2 in M 3 . Let Ric denote the Ricci tensor, 
let t denote the scalar curvature, and let L a b denote the second fundamental form. 
We let indices {i,j, k, 1} range from 1 to to and index a local orthonormal frame for 
TM; we let indices {a, 6, c} range from 1 to m — 1 and index a local orthonormal 
frame for TdM . On the boundary, e m will always denote the inward unit geodesic 
normal and ';' will denote the components of the covariant derivative. In Section^ 
we will perform a careful analysis of the local formulas involved and show: 

Lemma 1.4. There exist universal constants e 7 ^,^, so that: 

JdM 
JdM 1 
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T> /H,ai ,cx 2 \ V ,H I ^ c P/TC,a 1 ,Q 2 W )P / ^ c D/K,a 1 ,Q 2 ^ 

r ^,^,,^,,.^0 -ox-^^a /^o „o\__i3 
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T>/1Z,ai ,a 2 



O 1 ) +^a x ,a a W,P°> +4 7 ,a 1 ,a 2 <^ 1 ,P°) 

We omit the terms involving Sn for Dirichlet boundary conditions. 

We shall assume the following henceforth; it is properly part of Conjecture 11.21 
but we postponed it until we could introduce the appropriate notation. Note that 
it is crucial that we permit Oj to be complex so that the open set of Equation (|l.a[) 
is connected as this would fail if we only considered real oti . We refer to [2] where 
a related result was established when the specific heat is smooth: 

Conjecture 1.5. The junctions e^jm ai a2 appearing in Lemma \1.4\ have an ana- 
lytic extension to the connected open set defined by Equation U.a\) . 

1.6. The invariant (3q M . In Section [31 we will make a special case computation 
on the half-line to identify the constants £v/iz,ai,a 2 °f Lemma ll.4l Set 

f -1 if B = B v 
£v ' n '- \ +1 if B = B n 

Lemma 1.6. If (0:1,0:2) satisfies Equation \l.a\) . then 

- , , o-ai-q 2 „-l/2T ( 2-ai-a 2 \ r(l-q 1 )r(l-q 2 ) 



_ 2 _a,-a a7r -l/a r (tel) p( ai +a2 _ 1 y ^(l-ai) 



r(l-q 2 ) 

r( Ql ) 



The following result is an immediate consequence of Lemma 11.61 using the usual 
properties of the T function. We will give another proof in Section 2] that arises 
from certain functorial properties of the heat content asymptotics. 

Lemma 1.7. We have the recursion relations: 
_ 2(ai-2)( Ql -l) _ 

i-LJ £X>/TC,ai-2,a 2 — 5 £ T>/K,a 1 ,a 2 - 

6 — Oil ~ a 2 

_ 2(a a -2)(a 2 -l) _ 

\ z ) £ V/TZ,a l! a 2 -2 — 5 £ V/K,ai,a 2 - 

6 — Oi\ — Oi2 

/o\ c _ 2(a 1 -l)(a 2 -l) _ 

(,OJ tD/TJ. ai-l,a 2 -l — 5 £ K/V,a 1 ,a 2 - 

6 — Oil ~ Oi2 

Note that the roles of Neumann and Dirichlet boundary conditions are inter- 
changed in Assertion (3). 

1.7. Heat content asymptotics for Dirichlet boundary conditions. Dirich- 
let boundary conditions have been treated previously in [101 114] : we summarize 
those results as follows: 

Theorem 1.8. If (0.1, 0.2) satisfy Equation \l.a\l . then: 
0™ taa i<l>,P,DM = $ dM ev, ai , a2 {<t> ,P°)dy, 
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Pi^^i&PiDiBv) = $ dM {-j{ev, ai -i, a2 +e v , aua2 -i)Laa((f> ,P ) 

+Sv,a 1 -l,a 2 ((f> 1 ,p ) + £v, ai ,a 2 -l(^° , P^jdy, 
Pi^^aM^P.D.Bv) = J dM {~^{ev, ai -2,a 2 + £T), ai -l, a2 -l )L a a , P°) 

+£v, ai , a2 (E(j) ,p ) +£-d,q 1 -2,q 2 (0 2 ,/O°) + £v, ai M 2 -2{4> a ,p 2 ) 

,ai — l,a 2 — l~\~£T>,ai,a 2 —2j-Laa\ , P tP / 
+ (^j£V,a 1 -2,a 2 — \s~D,ui,a 2 -2 + \^T) ,a 2 ) {L a bL a b + Ric mm ) (</>° , p°) 

-£T>, ai ,a 2 {(i>% P%) +OT<j) p +ex), Ql _i :Q2 _i( ( ?!) 1 ,p 1 ) 

+ (^££1,0:1-2,02 + f £P,ai,Q 2 -2 + 'j£X',Qi-l,a2-l — j£T>,a 1 ,a 2 )L aa Li ) b(<t> , p°)}cfo/. 

1.8. Heat content asymptotics for Robin boundary conditions. The follow- 
ing is the main result of this paper; it will be established in Section @] using various 
functorial properties of these invariants: 

Theorem 1.9. If(a\,a2) satisfy Equation \l.a\l , then: 
C«»^aABb)= / e^ ai , a2 {^,P Q )dy. 

JdM 

>Y,V, „, (o.p. D. B K ) : / I e n , ai -i, aa (<l>\p )+e 1l , ai , a2 - 1 {<tf , ,p 1 ) 



JdM I 

_ Hsf^T £ ^^i- 1 ^2 + 57rx£K,oi,Q2-l}(-^aa0°, P°) 

+ {"SI 1 ^T £ K,a 1 -l,a2 - ■ Q ^T£K,a 1 ,a2-l}(5'TC'?!> ,p )|d2/ 

J a a/ I 



2 ,P°) +£TC, Ctl ,a 2 -2(^°,p 2 ) 



"2"{of^T £TC . ct i- 1 . ct 2-l + a2_2 e K,Qi,a 2 -2}(^aa0°, p 1 ) 



— — £n, ai ,a 2 {[R.ic mm +L ab L ab )(/) v , p v ) 

-. 2 1 „ 2 



2 3 — Ql- Q2 

1 1 aiQ2 , ,1/r r , , J,0 „0\ 



I r a i+°2- i I 1 am; -i/r r,,J.Q „0 

+ i 4 (3-ai-a 2 ) e7? -' Q i< Q 2 + 4 (ai-l)(a 2 -l) o rc,ai-l,a 2 -lM^aa^f>6 < / ) ,P 

-£ K , Ql ,a2«,P?a>+O-(T0 O ,p ) 

+e-R, ai .a 2 {E(jP,p ) + £ KjQ , 1 _i )Q!2 _ 1 (0 1 ,p 1 ) 

+ {(a 1 -l) 1 ( Q2 -l) £ K,«l-l,a2-l + 3-c t i-Q 2 e ^ ! "i^2}( S 'TC^ 'f ) 
+ {~^^2 £ rc,a 1 -2,a 2 ~ 7^^T e TC,Ql-l,a2-l}('S' , R,0 1 , P°) 
+ {~37^I £ K,a 1 -l,c<2-l _ ^^2£7Z,ai,a 2 -2}(»S'7?,0 , P 1 ) 

+ { 3-at-a2 £7Z < Ql < Q 2 + 5 (ai-l^aa-l) 6 ^"!- 1 ^-!} (^K^aa^ , P°) 

1.9. Outline of the paper. Our purpose is to a large extent expository so we 
shall give complete details rather than simply referring to previous papers. In 
Section [31 we discuss the local invariants we shall need. In Section [31 we make a 
special case computation on the half-line to determine (3q M . In Section [H we use 
functorial properties of these invariants to determine /3f M and /?f M - We must work 
in great generality in the context of operators of Laplace type in order to apply 
the functorial methods of Section [H It is one of the paradoxes in this subject that 
even if one is only interested in the scalar Laplacian in flat space, it is necessary to 
derive completely general formulas and then specialize them. 
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2. Local invariants 

By assumption, we may express the asymptotic coefficients Q2 as the inte- 
grals of local invariants over the boundary: 

P,D,B v/n )= [ P™ ua2 {^p,D,B v/n ){y)dy. (2.a) 

JdM 

Fix a system of local coordinates x — (x 1 , x m ). Set d x = (d Xl ) ai ...(d Xm ) am . We 
have 

P™ ua2 (<P,p,D,B v/K )(y) = ^ a2 {4> l ,p\dlg ih dlA\dlB,diSn) ■ (2.b) 

Here we need only consider a finite number of jets; we omit the Sn variables for 
Dirichlet boundary conditions and for Robin boundary conditions only differenti- 
ate S-jz tangentially. The local invariants are polynomial in the derivatives of the 
structures involved with coefficients that depend smoothly on the metric. 

2.f . A coordinate formulation. We summarize the arguments briefly - this also 
serves to motivate the power (1 +j — ai — a^jl of t in Coniecture lf .21 Fix a point 
P € dM and choose local coordinates y = (y 1 , y m_1 ) on the boundary which are 
centered at P. Let x :— (y,r) be the associated adapted coordinates. Let indices 
a, 6, ... range from f to m — 1. We then have 

ds 2 = g a bdy a o dy b + dr o dr where g ab = g a b(y, r) . 

We further normalize the choice of coordinates so that g a b{P) — Sab- This eliminates 
the smooth dependence on the metric tensor and ensures that the invariants of 
Equation (|2.bl) are polynomial in the jets of the structures. For a smooth function 
/ and a multi-index a — (ax, a m ) we let 

h = d a x \...d a xZ f. 

We define 

ordtey, s ) ■= \a\, ord(^) := \b\ + 1, ord(B c -) := \c\ + 2, 
ord(^) := £, ord(p e ) := I, ord(S , K (? ) := \d\ + 1 . 

Here, since SVj is only defined on the boundary, d only reflects tangential indices; 
we delete the g variables from consideration for Dirichlet boundary conditions 
as they play no role. 

Lemma 2.1. /3j^ a2 ((f>, p, D, B)(y) is homogeneous of weighted order j. 

Proof. Fix P £ dM and choose adapted coordinates so g a b{P) — Sab- Then the 
local formula for fif^ Q2 is polynomial in the variables 

{9ij/s, A^,B S , cj) e , p e , S K j}\ S \ >0 . 

We use dimensional analysis. Let c > 0, and consider the operator D c := c~ 2 D. It 
is then clear in a purely formal sense that e~ t<yC D > = e~( c ^ D and thus if T c is 
the temperature distribution defined by D Cl then 

T c (x;t) =T(x;c- 2 t). (2.c) 

We justify this formal computation by verifying that the relations of Equation (|I.cj) 
are satisfied. Suppose that T c is defined by Equation (|2.cl) . Then: 

(d t + D c )T c {x;t) =c~ 2 {d t T + DT}(x-c- 2 t) = and limT(-; c~ 2 t) = <f>(-) . 

We must check the boundary conditions are satisfied by T c ; this is immediate 
with Dirichlet boundary conditions so we set Bt> c '■= Bt> and ignore the subscript. 
The situation concerning Robin boundary conditions requires more work. We use 
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Lemma 11.11 to sec that D c and D determine the same connection. However the 
normal rescales; e m (c 2 g) — c~ 1 e m . Set: 

Bk c ■= V c -i Cm + S-n,c where Sn,c ■= c" 1 ^ . (2.d) 

We complete the proof that T c — e~ tDc cf) by verifying T c satisfies the rescaled Robin 
boundary conditions: 

B n T c = c- 1 {{V em +Sn)T{-,cr 2 t)}\ dM =c- 1 B n T(-,c- 2t )=0. 

Since g c — c 2 g, dx c — c m dx. We verify that the interior invariants j3^f rescale 
properly: 

/3?(4> )P ,D e ) = (-l) n /n\ ■X% g {(D n c <t> : p)} (2.e) 
= (-lf/n ■c m l° ieg {{{c- 2 D) n < P,p)} = {-lY/n ■ c m - 2n l a Reg {{D n ^ p)} 
= c m - 2n ^{^p,D). 
We apply Conjecture II .21 and use Equation (|2.ep to expand: 

oo 

p{4>,p,D c ,Bv/n c ){t) ~ Y, tncm ~ 2n Pn(<l>,P,D) 

oo 

+ Y. t(1+1 ~ ai ~ a2)/2 P™u^P'D c ,B v/ n c ). (21) 

3=0 

We apply Equation (|2.c|l to see: 



p{<t>,p,D c ,BT,/n c ){t)= (T c (x;t),p(x))dx c = c m (T(x; C - 2 t),p(x))dx 
Jm jm 

oo 

= c m /3(0, p, D, B v/n )(c- 2 t) ~c m Y. ic-H)»flf (0, p, D) 

n=0 

+c m ^(c- 2 t)( 1+ ^« 1 - a2 )/ 2 ^ )a2 (0 )/ ;,£>,Bp /7J ). (2.g) 



We equate powers of t in the asymptotic expansions of Equation (|2.fj) and Equa- 
tion (12. gD to see: 



^ iia2 (ip,i) c ,B B/Rc ) = c^+^-ipfX^^D^v/n) ■ (2-h) 
We use Equation (|2.a[) . Since dy c = c m_1 dt/, Equation (|2.h[) implies: 



P'DcBv/nJ = c'"- 1 / p?™ ua2 {4>,p,D c ,B v/nc ){y)dy 

JdM 

— c Pj,a 1 ,a 2 \ ( PiP-' 1J -' T>/1Z) 
„m-l-j+a 1 +a 2 I odM 



= c m-j- J+ai+a2 , ^ ia2 (0,p,D,B v/n )(y)dy. 

JdM 

Consequently (ignoring divergence terms), we have the following relation for the 
local invariants: 

/?|£U (4>, P, D c , B v/n< ){y) = C ^' +Ql+Q2 /3|f liQ2 (0, P, D, B v/n )(y) . (2.i) 

We wish to apply the formalism of Equation (|2.bp . We may expand: 

~ r- ai (0° + r c <pl + ...) - r- Ql (0° + r^ 1 + ...) . 

Equating terms in the asymptotic series, using the fact that r c — c • r, and arguing 
similarly with p yields: 

4>i = ^~Y and ^ = c Q2 -y- (2.j) 
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We renormalize the coordinate system and set x c — c ■ x; d Xci — c~ 1 d Xi . We 
introduce the coordinate systems x and x c into the notation as a computational 
aid; the local invariants a2 (y) do not, of course, depend on the choice of the 
coordinate system. We compute: 

g c ,ij(x c ,P) = c 2 gij (c~ 1 d Xl , c~ 1 d Xj )(x, P) = cr 2 c 2 gij{x, P) = . 

Consequently, g c ,ij(x c ) — gij(x) and we have that: 

9e,ij, aO^c; P) ^ ^ Qij .a(*£ ; P) • 

We expand: 

D c = c- 2 {g l ^d Xl d Xj +A k d Xk +B) 

= -(9 ij d Xc>i d Xe>j + c- l A k d XcM + c- 2 B) 
= ih ■ A% . ■ IU 

to see that A k (x c ,P) = cr 1 A k (x, P) and B c (x Cl P) = c- 2 B(x,P). Consequently 

A k j(x c ,P) =c- 1 -^Al(x,P) and B C;S (x c> P) = c' 2 -^ B s (x, P) . (2.k) 
Similarly, by Equation (|2.dp we have: 

Sn,cA^P)=c- 1 -^S 7Zd ix,P). (2.1) 

We use Equation Q2.jD , Equation f|2 .k[) . and Equation (I2.ip together with the 
observation that f3® M is bilinear in (0, p) to see: 

(h P> D c ,B v/n )(x c , P) (2.m) 

= /3f/ 1 , Q2 ( C al -V,c--V^ C - |S| 5^,c- 1 -l & Ul, C - 2 -l c 1B,, (2.n) 

^ A S nS )(x,P) 

= c-+-^ ja2 (c^/,c-V,c- |3| ^ )C - 1 - |S| 4,c- 2 -l c 1S ? , (2.o) 

c-^s^MP). 

We use Equation (|2.il) and Equation (|2.mp to conclude therefore: 

(c-V, c-V, c-" a l ftj7a> c^'^, cr 2 "^, c- 1 -!^-)^) 

= c-^2f liaa /A 4, B ft s Wi<r )(<?, P) • 

It now follows that flf^ a2 is homogeneous of weighted order j . There can be diver- 
gence terms which integrate to zero on the boundary and which are not controlled 
by this analysis. These terms play no role in the asymptotic coefficients and may 
therefore be ignored. □ 

2.2. A tensorial formulation. We have just discussed the homogeneity property 
of the invariants (3j from a coordinate point of view. We now consider the same 
property from a more invariant point of view. Fix a point P £ dM. Choose 
geodesic coordinates on the boundary centered at P. Form an adapted coordinate 
system for P in M. Then the only non-zero first derivatives of the metric are 
gab.m(P) — d Xm g a b(P)- The second fundamental form is given by setting: 

Lab ■— <?(V 'd*d y ,d r ) — T a bm = —\g a b,m ■ 

Keeping in mind the necessity to preserve the condition g a b(P) — S a b and observing 
that the Levi-Civita connection is unchanged by rescaling, we see that: 

Lab{g c ) = c 2 g{V c -i dya C~ 1 dy b ,c- 1 d r ) = cT Y L ah (g) . 
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Thus L is homogeneous of degree 1 under rescaling. Similarly we have: 

Rijki{g c ){P) = c 2 g((V c -i a[Ci V c -ia Xj - V c -i 9[Cj c -i d3:i )c~ 1 d Xk ,c~ 1 d Xl ) 
= c- 2 R l3kl (g){P), 

so that Rijki has order 2. Arbitrary partial derivatives of the metric at P may 
now be expressed in terms of these tensors and their covariant derivatives; there 
are universal curvature relations (see, for example, the discussion in [2H) but these 
play no role. Rather than considering tangential derivatives of S-r, we consider 
covariant derivatives of S-r with respect to the connection defined by D and the 
Levi-Civita connection on the boundary. Let $7^ be the curvature of the connection 
defined by D. Let ';' denote the components of multiple covariant differentiation 
with respect to the Levi-Civita connection of M and the connection V defined by 
D and let ':' defined similarly using the Levi-Civita connection of dM and V; we 
can only differentiate L and S-r tangentially. Each covariant derivative adds 1 to 
the order. The derivatives of the total symbol of D may then be expressed in terms 
of these variables so we may regard /3j jaija2 ((j), p, D,B)(y) as a polynomial in the 
variables 

{(f) , p , L aia2 - a3 _ ak , Ri 1 i 2 i 3 i i -.i 5 ...i k i E-.i 1 ...i k i Sfi:b 1 ...b k i ^ii i 2 ;«3 ■ ■ .** } 

which is invariantly defined and which is homogeneous of total order j. H. Weyl's 
theorem [32] then lets us write this in terms of contractions of indices; the structure 
group is the orthogonal group 0(m — 1) so the index m plays a distinguished role. 

2.3. The proof of Lemma 11.41 Lemma 11.41 now follows from Lemma 12.11 and 
the discussion in Section 12.21 once a suitable basis of Weyl invariants is written 
down. We integrate by parts and ignore divergence terms to replace the invariants 
JdM^aa, P°)dy and j dM {4P, p° aa )dy by - J dM (<f>° al p° a )dy. □ 

3. The proof of Lemma fOl 

We must determine the coefficient £r>/7j., Ql;Q2 of Lemma 11.41 which describes 
cf a 2 P> D, B). We begin with a computation on the strip: 

S := {(ati, 0:2) £ I 2 : ai < 1, CH2 < 1, o;i + a2 > 1} . 



The Dirichlet setting was discussed in [TOl [TT] so we concentrate on Robin boundary 
conditions; we may set S-r = as this plays no role in (3® M and simplify matters by 
considering Neumann boundary conditions. Let M — [0, 00), let D = —d x , and let 
Kj)/u{ x i V\ t) be the Dirichlet or Neumann heat kernel on M for t > 0. Let % be a 
plateau function which is identically equal to 1 near x = and which has compact 
support in [0, e) for e small. Set 

(f> ai ( x ) : = x~ ai x(x) and Pa 2 (x) := x~ a2 x{x) ■ (3. a) 

We shall use this notation subsequently as well. We have: 

1 if I = 1 , £ _ ( 1 if i = 
if£>0J an p ~\0 if £ > 

The heat content on [0, 00) is: 

/>00 f-OQ 

0(<l> ai ,p a2 ,D,B v /M)(t) = / K BT>/M (x 1 ,x 2 ;t)(j) ai (x 1 )p a2 {x2)dx 1 dx2- 

Jo Jo 

Then /3(cj) ai , p a2 , D,Bj\f)(t) is finite for all t > if and only if a\ < 1 and a 2 < 1. 
We also suppose that a\ + a% > 1. In [TUj we have shown that 

P(4> ai , Pa2 , D, B v ){t) = e v ^ a2 t (1 - ai - a2)/2 + 0(1) . 
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In the special case of the half-line we have that 

K Bn ( Xl ,x 2 ;t) = K Bd (x u x 2 ; t) + 2{^y 1 / 2 e-^ +x ^ . 
We compute: 

,p a2 ,D,Bjv)(t) = f3(4> ai ,p a2 ,D,B v )(t) 
+2{Ant)-^ 2 / / e~^ +x ^/ { ^x- ai x- a2 d Xl dx 2 (3.b) 



JO 

oo poo 



+2(4 7 rt)- 1 / 2 / / e-^ +x ^ /^x^ ai x-^{ X {x 1 ) X {x2)-l)dx 1 dx 2 . 
Jo Jo 

The second term in the right hand side above is evaluated by a change of variable 
x 2 = x\o '. By Tonelli's Theorem we have that it equals 

/>oo />oo 

2(4 7 rf)- 1 / 2 / / e-^ 1+x ^ 2 ^ t x^ ai X2° t2 dx 1 dx 2 
Jo Jo 

noo 

noo 
o- a -x\- a ^e-^ 1+ ^l^d Xl da 

= 2 l-a 1 -a a7r -l/2 r ( %~ai~a 2 ^ t{1 _ ai _ aa)/3 J°° ^ a2 ^ + ^ax+oa-a^ 

Note that if a± < 1 and a 2 < 1 then 

/•OO />1 

/ a- a2 {l + <r) ai+a2 - 2 dcr = / (cr"" 1 + o--" 2 )(l + cr) Ql+Q2 ~ 2 dcr . 
Jo Jo 



We will show that the third term in the right hand side of Equation (j3.b[) vanishes 
up to all orders. First note that there exists e > so x(x) = 1, < x < e. Secondly 
note that there exists a constant C < oo so < C, x > 0. Hence 

|x(*i)x(s 2 ) - 1| < (C 2 + l)l [e ,oo}(a!i)l[e,oD)(afa) • 

It follows that the third term in the right hand side of Equation (|3 .b[) is bounded 
in absolute value by 

/•oo />oo 

2(C 2 + 1)(4^) _1/2 / / e- l - x ^ x ^'^x~ ai x^ a2 dx 1 dx 2 



OO /"OO 



< (C 2 + l)t- 1 / 2 e- £ /< 4 *) / / e-W+^)/( 8 *'xr Ql ^ Q2 d a;i d a;2 
Jo Jo 

= 0(e -e 2 /(5t)). 

By LemmaO we have $% uQa {<l>,p,D,B) = [ s v/n , ai ^(4> a , p a )dy. The 

JdM 

above discussion permits us to decompose: 

^<B-D/rc' a i> a 2 — ^i3-D/Tz^ a i^ a 2 ~t~^ai,a2 wheie 

,2 — ot\ — a 2 . f 1 



2 -«i-«2 7r -i/a r (2 — ai — 02^ /■ ^_ ai +ff -« a ^ 1+(r )ai+a a -a dt7 
2 Jo 



— 0-Ql-a 2 ^- 



2 -«l-«2 7r -l/a r (2 Oi r ai +CT - Q2 )( 1 _ (T )a 1 +a 2 -2 dCT _ 
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The integral defining n ai ,a 2 clearly is well defined for ct\ < 1 and 0,2 < 1. By 
replacing cr by ^, we see that the integral with respect to a over [0,1] defining 
K aii „ 2 is equal to the integral with respect to a over [1, oo). Thus the integral with 
respect to a over [0, 1] is \ the integral over [0,oo). We may now use Gradshteyn 
and Ryzhik [35] (see 3.251.11) and Gradshteyn and Ryzhik [25] (see 8.3801 and 
8.384.1) to evaluate these integrals and establish Lemma 11.61 on the strip. By 
Conjecture 11.51 we may then use analytic continuation to establish Lemma 11.61 for 
the full parameter range given in Equation (jl.a[> . □ 



4. The proof of Theorem 11.91 functorial properties 

In Section 01 we determine the unknown coefficients a Q2 of Lemma 11.41 to 
complete the proof of Theorem 11.91 

4.1. Direct sum and product formulas. 

Lemma 4.1. The coefficients of Lemma \l.J\ are independent of the fiber dimension 
of the underlying vector bundle. 

Proof. Let V = V\ © V2 be a direct sum vector bundle over M, let D = D\ © D2, 
let = © 02, and let p = pi @ P2- With Robin boundary conditions, we take 
S-r = Stz.i © Sk,2- The problem decouples and consequently 

/3(0, p, D, B)(t) = (3{cj>x, p x , Dx, B){t) + P(<h,p2, D 2l B)(t) . 

This implies that: 

& ,<* (6 P. A B) = (3fX >a2 (0x ,px,Dx,B) + fl£ , Q2 (0 2 , P 2,D 2 ,B). 

The desired result now follows; again, we ignore divergence terms when passing to 
the local formulas. □ 

Suppose that (Mx,gx) is a closed Riemannian manifold of dimension mi and 
that (M 2 , 172) is a Riemannian manifold with boundary of dimension ni2- We form 
(M, g) := (Mi x M2,gi g 2 )- If Di are operators of Laplace type on Mi, we form 
D = Dx ® 1 + 1 <g> D 2 on V := Vx ® V 2 . 

Lemma 4.2. Adopt the notation established above. We have: 

(1) /3(0, p, D, B v/n ){t) = p{4>x , px , Dx)(t) ■ /3(0 2 , P 2,D 2 , B v/n )(t) . 

(2) fl£, a M>P> D >Bv/KK*i>V*) 

= E (-l) < 2lWl.Pl)(»l)-^a 1 ,a a (^.Pa,I33,iB I> / W )(« a ). 
k+2i=j 

(3) The coefficients of Lemma \1.4\ are independent of the dimension m. 

(4) £ T>/TZ, ai ,a 2 = £ ~D/TZ,a u a 2 , £ V/TZ,a 1 ,a 2 = ®> £ V/Tl,a u a 2 ~ ~ e V/K,ai ,a 2 ■ 

Proof. One has purely formally that 

e ~ tD ^/n = e -*Di fa g, e -iD 2 , Bv/ll ^ ^ 

We verify Equation (|4.a|) using Equation (|l.c[) as follows. Let Ti be the temperature 
distributions defined by 0j on Mj. We compute: 



(fit + D)(Tx ® T 2 ) ={(d t + Dx)Tx} ®T 2 +Tx® {{d t + D 2 )T 2 } = 0, 
Iim{Ti® T 2 )(-;t) = lim{Ti(-;t)} ® lim{T 2 (•;*)} = 0i <g> 2 

#D/tc(01 ® </>2 ) = 01 <8> {S-D/7J.02} = . 
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Assertion (1) is now immediate. Since the boundary conditions play no role on the 
closed manifold Mi, we have 

W.p.^W-Vf^- / (D n 4>i, pi)dxi . 

Assertion (2) now follows from Assertion (1) by equating terms in the asymptotic 
expansions. 

We argue as follows to prove Assertion (3). Let Mi = S 1 , let D± = -d%, and let 
4>i — pi = 1- Since T±(xt;t) = 1, 

P{4>u Pl ,D 1 )(t) = 2ir. 

We use Assertion (2) to see that: 

f3™, a2 (^p,D,B v/n ) = 2n-pfXt a2 (^,P2,D 2 ,B v/K ). 

This implies that for the local formulas it does not matter if we compute on M or 
on M2. Consequently the coefficients of Lemma |4~21 are independent of m. 

We complete the proof by considering Assertion (4). With Robin boundary 
conditions, we take S-jz = as it plays no role in the evaluation of the coefficients 
we are considering. Let (Mi,<?i,<^i,pi,-Di) be arbitrary. Let M2 — [0,7r] with the 
usual flat metric, let D2 = and let <j> ai and p a2 be as given in Equation (|3.a[) . 
Since the structures on M 2 are flat, 

odhh (A „ n n \ — j £ V/n, ai ,a 2 if j = 
Pj, ai .a 2 { ( P2,P2,-L>2,B T , /R ) - < Q ifj>0 

Consequently Assertion (2) yields: 

/^^(^P- A#d/tc) = -£v/n. ai ,a 2 / (Dtfa, pi)dxi 

J Mi 

= £n, ai . a2 (4>l- t aa +0-Tl(j)i,pi)dx 1 (4.b) 

J Mi 

= £v/n., ai .a 2 J |- {<f>l;a, Pl;a) + {E(/>1, pi) + (OT01, pi)| dxi . 



On the other hand, when we use the formulas of Lemma 11.41 directly, most of the 
terms vanish and we obtain: 

P. A BtO (4.c) 

= / |4 2 /R.oi,02^i«'/'l;«) + (4>/W,ai,a a - B ^l+ e I>%,a 1 ,a a r ^l»/ , l)} <fa l ■ 
J Mi k ' 

The desired result follows by comparing Equation (|4.b[) with Equation (|4.c|) . □ 
4.2. Index shifting. The following observation is simply change of notation: 
Lemma 4.3. 

(1) Suppose that <p° = 0. Then cf> can also be regarded as being an element 
4> € IC ai -i(V) where = <j) i+1 . We have: 

Pf™i, ai -i, a2 {I P, D, B v/n ) = fi% uaa {<j>, p, D, B VJn ) . 

(2) Suppose that p° — 0. Then p can also be regarded as being an element 
p £ JC a2 -i(V*) where p l = p l+1 . We have: 

fl?i, auaa -M P, D, B v/n ) = fi% uaa (</>, p, D, B VJn ) . 
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(3) We have 



t X>/K,ai 


,OC2 


— £ D/K,qi- 


-1,Q2 ' 


f 3 

c I5/7?.,ai,a2 


— £V/K,ai 


,«2 


-1j 


t B/R,ai 


i<*2 


— £v/Tl,at- 


-2,CK2 5 


F 7 


= £X>/TC,ai 


,«2 


-2> 


5 

t X>/TC,ai 


|C*2 


- r 2 


-1.Q2 ' 


f 8 

c I3/7?.,ai,a2 


- f 2 


,CK2 


-1' 


14 

t X>/TC,ai 


|C*2 


— £T>/Tl,ai- 


-1,012 — 1 7 


7?.,ai,a2 


p 15 


) 




F 18 




-15 
















t K,ai,Q 2 -l- 













Proof. Suppose 0° = 0. Then 

0^ r - ai (r^ +r 2 2 + ...) - r - ^ 1-1 ^ 1 + r0 2 + ...). . 
Consequently, <j> € /C ai _i. We expand: 
p(4>,p,D,B v/n )(t) 

oo oo 

= p(4>,p,D,Bv/n){t) 

OO oo 

n=0 k=0 

We set fe + 1 = j and equate terms in the boundary asymptotic expansions to prove 
Assertion (1). The proof of Assertion (2) is similar and is therefore omitted. 

We apply Assertion (1) to prove Assertion (3). Suppose <f>° = 0. By Lemma [Ql 

JdM 

C^A/'-AM 3 / cv/K, ai -i, a2 (4> Q ,P°)dy. 

JdM 

Consequently, Sj-,, n Qi Qa = £vfR.,ai-i,a 2 - The remaining assertions of the Lemma 
are established similarly. □ 

4.3. Relating Robin and Dirichlet Boundary Conditions. We work on the 
interval [0, n]. Let ^ c e R. Set 

A:=d x + c, A*:=-d x + c, D = A* A = AA* = -(<9 2 - c 2 ) . 

Let £?x> be the Dirichlet boundary operator and let Bn<t> = define the boundary 
condition A*4>\qm = 0; this is the associated Robin boundary operator: 

B v ^ := 0(0) © 0(tt), B n <\> := (<£'(0) - c0(O)) 8 (-0'(tt) + c^tt)), 
S n (0) = -c, S n {ir) = + Cl E = -c 2 . 

Lemma 4.4. If$t(ai) < — 1 <™d 3?(a2) < — 1, fftera: 

(1) a t /?(0, p, D, = -0(D^ p, £>, Bp/ W )(t). 

(2) d t ^,p,D,B v/n ){t) = -/3(<fi,Dp,D,B v/n )(t)- 

(3) W,p,A^)(*) = -^(AV,A>,£>,Bp)(i). 

(4) 9 t /?(0, p, D, B v ){t) = -p{A<t>, Ap, D, B n ){t). 

(5) ^,Z+2,a 2 (D<P,P,D,B v/R ) = - - ai - a 3 )0f£ uaa (<t>,p,D,B v/1i ). 

(6) Pl%, a2+2 (<t>,Dp>D,B v/n ) = -i(l + j - ai - a 2 )/8^, a2 (AaABd/r). 

(7) $£ + i, aa +i(^ ^*ft A Bo) = -K 1 + J - «i - <X2)Pf£ >aa te, P, D, B n ). 
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(8) f3fX+i, a2+ i(^MD,B n ) = -U 1 +j-^-^)PfX, a MP^ D ^^- 

(9) £ TZ,a u a 2 ~ t^ai-l^ 11 !- 1 .^ _ Q2 - 1 £ K,Qi ,a 2 - 1 } • 



(io) 4V 



0:2 



ai- 


1 


1 






2 


1 





( 12 ) £ K, 0l ,Q 2 — (ai-lKas-l)^" 1 " 1 '" 2 " 1 ^ 3-ai-a 2 eK ' ai ' Q2 

Proof. We begin by establishing the spectral resolutions of the Dirichlet realization 
D-p and of the Robin realization D-jz. For n — 1,2,..., set: 



4>n '■— y f sin(na;) and A„ := n 2 + c 2 . 

Because {4%} is a complete orthonormal basis for L 2 , because D<ff^ = \ n <$^, and 
because <^(0) = (j^i 71 ) = 0, we may conclude that {<^,A„}„ e N is the Dirichlet 
spectral resolution of D. Similarly, set: 

{(fi^} is a complete orthonormal basis for L 2 with D<f% — X n <P^- Since 

A*<ft\ BM = \- n x l 2 A*A<fc\ dM = X-^DcfiloM = \l/ 2 <f>v\ dM = 0, 

B-jzA^ — so the eigenfunctions satisfy Robin boundary conditions. Conse- 
quently {(ffi, A„}„ e N is a Robin spectral resolution of D. If <j> e i 1 , let 

:= [" <Kx)<%'*{x)dx 
Jo 

be the associated Fourier coefficients. We then have 

oo 

(3(<P,p,D,B v/n )(t) = J2 e ~ tXn ^ /U (^ /1Z (p), 

n=l 

oo 



n=l 

Since 3J(ai) < —1 and 5i(a2) < — 1, we may integrate by parts to see 

V/Kf 

— i 

M JM 



In 



1 {D4>) = [ Dcj) ■ <f%/ n dx = [</>■ Dfl'Kdx 

JM JM 
JM 

Assertions (1) and (2) now follow. To prove Assertion (3), we note: 
^(A*<f>) = / A**- <%dx = /V Atfdx 

JM JM 

= / <f> ■ (j>n<lx = ^Kjn(<P) ■ 

JM 

Thus we may establish Assertion (3) by computing: 

oo 

d t /3(^,p,D,Bn)(t) = £-A„e- tA " 7 ^) ■ 7?(p) 

n=l 

oo 

= -E z~ tXn ln{A*ct>) ■ jV { a* p ) = -0{A*<f>, A*p, D, B v )(t) . 

n=l 

Assertion (4) follows similarly once we observe that {A* <j>% / y/X^ , A„}„ e z is du- 
ally a spectral resolution of the Dirichlet Laplacian; the roles of A and A* being 
interchanged. 
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We now establish Assertions (5)- (8). Since D(f> € K. ai+ 2(V), we have: 

oo 

d t p(<t>,p,D,B v/K )(t) ~ £n- (-lrt^/nlX^UD^p)} 

n=0 

1 °° 

+ 2 EC? + 1 - «i - «2)t (j - 1 - Ql - Q2)/2 /3|^ 1 , a2 (0, ft A Sp/r) 

OO 

= -(3(D0,p,D,B v/n )(t) ~ ~^(-l)V/«I Rcg {(^+V,p)} 

oo 

_ £ #+Hai+2) _ a2)/ ^M +2)aj p) ^ ^ 

fc=0 

The term with n — plays no role. Setting n = £+1 equates the summation involv- 
ing the regularized integrals. Setting k = j and equating terms in the asymptotic 
series for the boundary terms establishes Assertion (5); the proof of Assertion (6) 
is similar. To prove Assertions (7) and (8), we may integrate by parts to express 

l Rcg {(D n A*cj ) ,A*p)}=l Rcg {(AD n A*ct>,p)}=l Rcg {(D n AA*ct>,p)} 
= l Kes {{D n+1 cf>,p}}, 

l Rcg {(D n A^Ap)}=l Rcg {(A*D n A^p)}=l Rcg {(D n A*A0,p)} 
= l Rcg {(D n+1 ^p)}. 

The remainder of the argument now follows similarly. 

To prove Assertion (9), we assume 5J(ai) < — 1 and 5R(a2) < — 1. We then 
use analytic continuation to obtain the result on the full parameter range. Adopt 
the notation of Equation (I3.ap and set <f> — 4> ai , and p = p a2 . Only the value 
at r — is relevant; integration over the boundary is just evaluation at 0. Since 
D = — (<9 2 — c 2 ), the metric and associated connection are flat. We have: 

A*(j> = {-d x + c){r- ai )=r- a ^- 1 {a 1 +cr}, 

A*p = {-d x + c)(r- a2 ) = r- a2 ~ 1 {a 2 + cr}, 

<f>° = l, </> 1 =0, p° = l, p 1 =0, 

(A*<j>)°=ax, {A*<t>Y=c, (A*p)° = a 2l (A*p) 1 =c, 

<t> 2 =0, p 2 = 0, (A*p) 2 = 0, (A*^) 2 = 0, 

Sr(0) = -c, E= -c 2 . 

We use Assertion (7), Lemma fl~4l and Lemma l4~3l (3) to see: 

fl% 1 +i,c*+i(A*<l>,A*p,D,B 7 >) 

= ^, ai+ x^+x{A*4>) 1 ■ (A*p)° + 4,a 1+ i,a a+ i(^V)° • [A*pf 

= C ■ {a2£x>,ai,c*2 + l + Oti£T>,a\+l,a^} 

= -1(2 - a x - a 2 )/3f A Q / liQ2 (0, p, D, B n ) 
= |(2 - ax - a 2 ) ■ ce^ >aua2 . 
We solve these equations to establish see: 

£ K,a u a 2 = 2-ai-a 2 ( a 2 £ X>,ai,a 2 + l + <*l£'D,m+l,ai) ■ 

Applying Lemma [TTTl (3) then yields: 

..15 _ r i F l _ i 

t R,Qi,a 2 ~~ l a 1 -l o7 ^ ! a i- 1 ^ Q 2 ct2 _i t 'R.,Qi,Q2-lJ 

Assertion (10) and Assertion (11) follow from Assertion (9) by using Lemma R~31 to 
index shift. We take j — 2 and examine the coefficient of c 2 to complete the proof 
of Assertion (12) by computing: 
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PiM +1 , a2+ i(A*4>,A*p,D,B v ) 

= (A*(j)) 1 (A*p) 1 s]^ ai+l a2+1 - c 2 (^*0) o (-4*p)°e|) :Ql+ i iCe2+ i 

2 2 

— C • £D J Q 1) a 2 ~~ c a l Q; 2£T),ai + l,a2 + l 

= _i( 3 _ ai - a 2 )p$M ua2 (<t>, p, D, B n ) 
= -i(3 -a x - a 2 )c 2 {e^ Qi Q2 - e TC , ai , Q J- 
We solve these equations and use Lemma 11.71 to see: 



.16 2 , 2a 1 a 2 _ , _ 

'•R,ai,a 2 — 3-ai-a 2 X) ' Ql ' Q2 3- Qi -q 2 C ' Ql + 1 ' Q2 + 1 ^ t K,Qi,a 2 

2 3-ai-a2 

~~ 3-ai-a 2 ' -2(ai-l)(a 2 -l) ^.ai— l,a 3 — 1 

2aia 2 l-ai-a2 f - _i_ 3— ai — a 2 - 

1 "3~Qi-a 2 -2aia 2 fc, R-,ai,a 2 T 3_ ai _ Q2 =K,ai ,a 2 

= (ai-lKaa-l) 67 ^."!- 1 ^ 2 -! + 3_ ai _ a2 £K,ai,a 2 ■ 



a 



4.4. The proof of Lemma 11.71 Lemma \A. 31 (3) focuses attention on the shifted 
invariants £s, ai — 2,a 2 > £ B,a 1 ,a 2 -2, and £b, Qi -i, Q2 -i- Lemma ll.7l expresses these 
invariants in terms of the fundamental invariants £B, ai , a2 . These relations follow 
from the standard properties of the T function. It is instructive, however, to use 
Lemma 14.41 to establish these properties as an independent check on our work. 

We adopt the notation of Lemma 14.31 with c = 0. We use Equation (13. a[) to 
define <fr ai and p Q2 on [0, n]. We take j = and apply Lemma H~4l and Lemma l4~4l 
to see: 

= ^/ ? 0,fi f t +2,Q 2 (o ! l(Q ! l + 1)051+2, Pa 2 i D^v/tz) = -Si (Si + l)£xyR,6i+2,a a 
1 — 5i— a 2 a3M /jt n R , ^ 1 — «i— a 2 

Consequently: 

— Si(Q!i + l)eX)/K : ai+2,a 2 = ~ 1 "2 ° 2 £ P/R,tti,a 2 ■ 

Lemma [T771 (1) now follows after replacing setting <Si — a.\ — 2 and &2 — Q.2- The 
proof of Lemma ll~7l (2) is similar. To prove Lemma [1771 (3), we set c = to obtain 
A = —A*; the roles of T> and TZ are then symmetric in Lemma 14.41 so we may 
compute: 

= aiQ;2/3o,i f 1 + l i 5 1 + l(0Q 1 +l,PQ 2 +l,Sl7/TC) = "l"2£x>/TC,ai+l,Q2+l 

_ l-ai-a; fl9M / 'j. _ n fc? \ _ 1— fir- ag 

— 2 P0,5i,a 2 l < ? ) >P^> O TC/-Dj — 2 £ TZ/V,a 1 ,a 2 

Lemma fTTTl (3) now follows after setting Sj = a; — 1. □ 



4.5. Warped products. Let T m_1 = [0, 27r] m ~ 1 denote the torus where we iden- 
tify ~ 2tt to form a closed manifold. Let {y\, y m -\) be the usual periodic 
parameters where < yi < 2tt, and let 

A/ := T"'- 1 x [0, tt] . 

Let / a € C°°([0, 7r]) be a collection of smooth functions which have compact support 
near r = with / a (0) = 0. Let xi r ) G C°°([0,7r]) be the mesa function described 
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previously; \ is identically 1 near r = and identically near r — 1. Set 
ds M = X e2fa(r) dya ° rfy Q + <ir o dr, (f> 2 := y(r)r _Q;i , 

a 

£>M := - 5^ e_2/a(r) C - P2 == X(r)r-«= , 

a 

D 2 :=-d%, <p M (y,r) := <p 2 (r), PAi(y,r) = e~~£« fa ^ p 2 {r). 

Let 2? = B-d or let £?k = d r + S-jz,o- 

Lemma 4.5. Adopt the notation established above. Then: 

(3 3 , aua2 ((f>M,PM,D,B)(y,0) = f3^ auOl2 (<j) 2 ,p 2 ,D 2 ,B)(0) . 
Proof. Let T 2 := e- tD2 - B (f> 2 and let T M := e- tDM - B 4> M . We show that 

TNi{y,r;t) = T 2 (r;t) 
by verifying that T satisfies the defining relations of Equation (| 1 .c[) : 

(d t + D M )T 2 = (dt + D 2 )T 2 = 0, UmT 2 (-;t) = 2 =0 M , BT 2 (-;t)=0. 

Since the volume element on M is given by e^»^ a ' r 'dydr, we complete the proof 
by computing: 

P{cf> M ,pM,D M ,B)(t) 

= /t»-ix[o,i] T2(r; <){e_ Eo /o(r) ^ 2 (r)}eEa /a(r)d y*- 

= /T^x[o,i] T 2( r ;^W¥^(2^r7 [ o,i] r 2(nt)p2(r)dr 

= (2ir) m - l /3(<t> 2 ,p 2 ,D 2 ,B)(t). □ 

4.6. The proof of Theorem ll.9l We have computed most of the unknown coeffi- 
cients in Lemma fO)! We use adopt the notation of Section FOl and apply Lemma l4~51 
to compute the remaining coefficients; the recursion relations of Lemma If .71 also 
play a crucial role as do our previous computations. We adopt the notation of Equa- 
tion (13.a[) and consider the following example. Let <p 2 := r~ ai x and p 2 :— r~ a2 \ 
where \ is identically 1 near r — and vanishes for r > |. We then have: 

4 = P \ = 1, 4 = pi = for I > 1 . 

Consequently there exist suitably chosen constants {cj} (which play no role in our 
development) so that: 

Pj,a 1 ,ot 2 (4>M,PM,D M ,B)(y,Q) = Pj iai!a2 (<f> 2 ,p 2 ,D 2 ,B)(Q) 

= c i ■ s n,o for > ■ 
We shall ignore this term and concentrate on other terms. Let uj be the connection 
f-form of the connection on V defined by D and let uj be the dual connection 1-form 
of the dual connection on V*; uj + uj = 0. We let e a '■= d Va for 1 < a < m and 
e-m '■— dr] this is an orthonormal frame on dM. We use Lemma If .f I to compute: 

r h — -fs h e 2 [ f» r i, m — -fp 2 f a 8 h 

i abm — J a u ab c 5 ab — Ja c u ab: 

i amb — J a u ab^ i 1 am — J a u a> 

L a b \dM = r a 6 m |aJVf = — fa$ab\g M > 

u a = 0, uj a = 0, 

U m = ~k fa) w ™ = = 2 /a • 

Consequently: 

Rambm\@M — 9{{^a^7n V m V fl )efc, 6 m )|gj^ = {T ac -f mb ab }|<9M 



HEAT CONTENT ASYMPTOTICS 



19 



\3M 



RiC mm | aAf — J2a {fa + (fa) 2 }\g M i 
E\dM = { — d m UO m — Ojf n + LO m Faa m } 

~~ "[^Sa/a — 4 Ea,6 f af b 2 Sa,i) /a/b 1 
= {lEa^o + I E a ,b /a/&} 



At this stage, the connection 1-form enters in a crucial fashion as we must use the 
dual connection and the dual connection 1-form Co to covariantly differentiate p. 
Thus Vg r V = (d r ip + w m V) if ip e C°°(V) while V ' dr ip = (d r ip - u m tp) if ip € V*. 
We use Equation (jTTbJ) to compute: 



{V 9r (r^ M )}| 9M = {(^-5Ea/i)(l)}U=-|E„/i| a ^ 
H(V a J 2 (r a ^M)}|aM = \{{dr-\Y, a fa?^)}\ dM 

8 "[X/a.h / a/ b ~ 4 Ea /a 1 



9Af 



Pm 
Pm 



\i)M 



= {V* dr (PM)}\ 9M = {(dr + \ Ea fa){e- ^ f «)}\ dM = -{hY Ja fa} 

|{(vs,) 2 pm} 



1 3M ' 



dM 



H(^ + |E a /^) 2 (e-^^)}| 



9M 



8 |Ea,i)/a/ii 4 So^n} 



AM 



To ensure that the Robin boundary on M takes the form B = d r + Sn^o, we set 

#K = -5^,0 - LJ m \ dM = S^q + { 2" Ea /a}|aM ■ 

4.6.1. TTie coefficient ofJ2 a fa in P\ M ■ 



0~ 2 £ TC,ai,a 2 e R,cti,Q 2 2 £ TC,ai,a 2 + 2 e; R.,al,a 2 ■ 



l c •3 



1_15 



Consequently 



" 2 fc 7?.,ai,ct2 



I =-3 



_i_ 1 P 15 



2 c R,ai,aj 2 7J,ai,a 2 



2^ ^71,01 — 1,02 ^7Z.,ai,a2 — 1 { a , _l ^71,ai — l,et2 a 2 — 1 ^7?.,ai i Q 2 — l}} 

1 f ttl - I t»2 c 1 

— 2lai-l t K,ai-l,ct2 "T" a2 _i t 'R.,ai,a 2 -l J" 

We use Lemma T4. 31 to shift indices and compute: 



F 5 = F 2 

1Z,oti 7 a2 7t,ai — l,at2 

— 1 f ai-1 , j a; 

2 I Ql _2 t K ! "i-2."2 T" Q2 _ 

,-8 _ ^2 

7Z,Qi,a 2 7?.,ai,a 2 — 1 



1 f eti I a 2 — 1 - 

~~ 2 lai- 1 t 7J.,ai — l,a 2 — 1 t Q2 _2 t TZ,a 1 ,a 2 -2 



I £lZ,a%— l,a2—l }) 
}• 



4.6.2. T/ie coefficient of Sn,o E a /a * n 02 



DM 



-16 _ 1 ,.17 _ 1 ,.18 _ ,.19 

'7J,Qi,a 2 2 R,ai,Q2 2 7i,ai,a 2 R,ai,Q 2 • 



Consequently 

-19 

7^,01,02 



-16 



(ai-l)(a 2 — 1) 



I«,17 



^TZ.Oii — l,a2 — 1 



1 _18 



ai-2 



-if 



_9 £-7^.,ai — 2,Q2 



3 — Q£l — a 2 
J3l £7£,a i — 1 ,a 2 — 1 } 
^Z2^,ai,Q2-2} 
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_ 1 2 f 9 , (ai-2)(g 1 -l) (g 2 -2) (g 2 - 1) ] 

— 2 3-qi-q 2 \ cti-2 T a 2 -2 J fc ^»ai.«2 

+ 2 1 ( ai -l) 1 (a 2 -l) { 2 + ("1 + ("2 ~ l^e^ai-l.aa-l 

_ ai+a 2 , 1 ai+«2 - 

~~ 3-Qi-Q 2 tTC . Q i. Q 2 ' 2 (Qi-l)(a 2 -l) tK . a i- 1 ' Q 2-l- 

4.6.3. T/ie coefficients ofj^af'd and J2a(fa) 2 in fil M ■ We obtain the relations 

fl--Ip 4 -4- Ip 6 - lr 7 - r 9 and 

u ~~ 4 fc TC,ai : Q 2 + 2 fc TC,ai,a 2 4 fc TC,ai,a 2 fc 7J.,ai,a 2 ' tM1U 



f) — — r 9 -4- F 11 



This shows 



P 9 = P 11 = — -F A A- -F 6 — -P 7 

t, 1Z,a 1 ,a 2 t, 1Z,a 1 ,a 2 4 1Z,a l ,a 2 ' 2 R,ai,aj 4 R,»i,a2 

4^7?., ai — 2,a 2 4 £-7^,ai ,a 2 — 2 ~i~ ^^TZ,ai,a 2 

= f 1 2(ai-2)(ai-l) _ 1 2Q 2 -2)(a 2 - 1) ill,. 

\ 4 3-ai-a 2 4 3-ai-a 2 > 2J fc7 ^ Q i> a 2 

_ 1 ct\— 2ai+a 2 — 2q 2 + 1 

— — 2 3-ai-a 2 e K : ai,Q 2 - 

4.6.4. The coefficient of^J'Ji m f3$ M . 

n — Ic 4 4. ic-S -4- Ic-6 1 I P 7 , 1 _8 1 _10 

u — 8 G R ,ai,a 2 ~ 2 7?-,Q:i,Q!2 4 "R-,ai,a: 2 ' 8 T^,Gti,ot2 2 *R*,cti,ot 2 ~ lZ : a\,cn 2 

,1 14 , 1..16 _ 1-.17 _ 1_18 _ l p 19 

~ r 4 R,ai,Q 2 4 u K,ai,a 2 4 Ti,ai.a 2 4 K,ai,ai 2 'R.,ai,a 2 ' 

This implies 

c io _ _i_ _ i_ _ 1. _ i_ 

t 7?,,ai,a 2 — g fc 7?.,ai— 2,a 2 4C7?,,ai,a 2 gtR iQl . a2 _2 4 ^TZ.ax — l,a 2 — 1 

2"! 2"L Q *_2 e "R.,ai-2,a 2 + a2 ^-i e 'R.,ai-l : Q 2 -l}} (~2 £ K,Qi,Q 2 ) 

— 2"t~ 2 { q 1 -1 £ 'R^i-1,q 2 -1 + ct 2_2 £ 'K,ai:Q2-2}} (~ 2 £ TC,ai ,a 2 ) 
~l{ (qi-1)(q 2 -1) £k ^i- 1 ^2-1 + 3 „ Ql _ Q2 £ K,ai,Q 2 } ( — 4 £ TC,ai ,a 2 ) 
+ 4"t~ Ql -2 £ ^ Q i-2> Q 2 — a2 _i e 7?.,ai— l,a 2 — l} (4^^01,02) 

+ 4"t~ Ql _i £ 7Z,ai-l,a 2 -l — a2 _2 £ K,Qi,a 2 -2} (l £ 7Z,ai,a 2 ) 

1 1 f ai+a 2 1 1 a 1 +a 2 \ /Ic!9 "1 

^2l3-t>i-a 2 ^ a i. a 2 f 2 (Qi-l)(a 2 -l) t7 ^^i- 1 ^2-l/ V 2 t TC,a 1 ,a 2 i 

_ f 1 I 1 1 1 1 c , r lil a 2 -l 1 1 1 

\ 8 T 4ai-2 4 ai-2 J t K,a 1 -2 : Q 2 "T ^ 8 t 4a 2 -2 4 a 2 -2 J t K,"l :"2 -2 



■{■ 



1 , 1 a 2 , 1 ai _ 1 1 _ 1 1 _ 1 1 

4 ' 4a 2 -l 4ai-l 4 (ai-l)(a 2 -l) 4a 2 -l 4ai-l 

I 1 ai+a 2 \ 
r 4 (ai-l)(a2-l) J t K,ai-l,a 2 -l 



f_ 1 _ 1 2 1 1 ai+a 2 ] 

1 4 4 3 — cci— a 2 2 3 — Qi— a 2 J 

ic- _4_ ic- 4- 3ai+3a 2 — 5 ,. 

g t R,Qi-2,a 2 T g^K,Q!i,Q 2 -2 "T" 4(3 — Ql — Q2 ) b 7l,ai,a 2 

\ {-1 + 1 + 1 + ( Q "ij)"a7-1) } E ^,ai-l,a 2 -l 



_ (ai-2)(ai-l) + (a 2 -2)(ai-l)+3ai+3a 2 -5 „ , 1 a\a 2 

~ 4(3-Qi-Q 2 ) t.n, ai ,a 2 -T 4 ( Ql _ 1 )( ct2 _ 1 )fcK,ai-l,a 2 -l 

_ nj + a 2 -l , 1 «m 2 

— 4(3-Qi-Q 2 ) t ' R '. a i^2 f 4 (Qi-l)(a 2 -l) t7 <. Q i- 1 . Q 2-l- 

We have determined all the unknown coefficients in Lemma 11.41 This completes 
the proof of Theorem 11.91 □ 
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